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Abstract 

We compute the decay constants for the heavy-hght pseudoscalar mesons in the 
quenched approximation and continuum hmit of lattice QCD. Within the Schrodinger 
Functional framework, we make use of the step scaling method, which has been pre- 
viously introduced in order to deal with the two scale problem represented by the 
coexistence of a light and a heavy quark. The continuum extrapolation gives us a 
value = 192(6)(4) MeV for the Bg meson decay constant and /^^ = 240(5)(5) 
MeV for the Dg meson. 

Key words: Heavy flavors; decay constants; 6-physics; lattice QCD 



1 Introduction 

The amount of CP violation occurring in the in the Standard Model de- 
pends upon the unitary Cabibbo-Kobayashi-Maskawa (CKM) matrix. This is 
explored through the so called Unitarity Triangle Analysis [1,2,3] that requires 
a deep theoretical understanding of the B mesons decay properties, given the 
increasing accuracy of the 5-factory experiments (for a recent review see [4]). 
A crucial quantity, the B meson decay constant Jb, has already been calcu- 
lated by means of different techniques: Heavy Quark Effective Theory (HQET) 
[5,6,7] QCD Sum Rules [8,9,10,11,12,13,14] and lattice QCD (LQCD). 

The latter moves from first principles and has to face the problem of prop- 
erly accounting two largely separated energy scales, i.e. the heavy and light 
quark masses. This imposes stringent limits on the values of the lattice bare 
parameters: the lattice spacing has to be small enough to allow a good descrip- 
tion of the highly localized heavy quark, and the number of lattice points has 



to be large enough to accommodate the widely spread light quark. A direct 
calculation would require lattice sizes of the order 100'^ points. These limits 
are too much demanding for present supercomputers, even in the quenched 
approximation, and in the case of full QCD also for the next generation su- 
percomputers. To overcome these difficulties, different strategies have been 
adopted (see [15,16] for recent reviews). 

One way consists in working with propagating heavy quarks with masses 
in the region of the physical charm quark and to extrapolate the results to 
the 6-quark, using HQET scaling laws [17,18,19,20,21,22]. The systematic er- 
rors in these calculations are dominated by the uncertainty on the functional 
form to be used in the extrapolation, depending in turn upon the limits of 
validity of HQET when applied to the c-quark. Another method comes from 
the non-perturbative matching of lattice QCD and HQET in the static limit 
[23,24]. Other possibilities are offered by the Non Relativistic approximation 
(NRQCD), where the heavy quark mass expansion is taken in the operators 
as well as in the Lagrangian [25,26], and by the Fermilab approach [27], where 
one expands in the mass of the propagating heavy quark in lattice units. All 
these methods give results with systematic errors that makes them fully com- 
patible within the error bars [15,16]. 

In a previous publication [28] we have introduced a new method to per- 
form the determination of /b and, more generally, to face two scales problems 
in lattice QCD. The so called step scaling method [15] has been applied to give 
a first numerical result for this quantity at fixed lattice spacing and, in [29], 
to perform the first calculation of the 6-quark mass in the continuum limit of 
quenched lattice QCD. The idea behind the method consists in using a QCD 
propagating 6-quark on a small volume, calculating the finite volume effects 
on the heavy-light decay constants and, finally, using the very mild depen- 
dence of these effects upon the heavy quark mass, to obtain a final result in a 
large volume. 

In this paper we extrapolate our results to the continuum by repeating the 
various steps of the previous calculation at different lattice spacings and fixed 
physical quantities. The major assumption undergoing our method, i.e. that 

the finite size effects on an heavy-light observable have a milder dependence 
upon the heavy quark mass than the quantity itself, is shown to hold in the 
continuum limit. 

The general features of the method are described in sec. 2, in sec. 3 we set 
the notations, the numerical results are presented in sec. 4. Conclusions are 
drawn in sec. 6. 
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2 Overview of the method 



A detailed explanation of the method can be found in [29]; here we shortly 
review the basic features in order to set the notations. 

A non-perturbative determination of the heavy-light decay constants from 
lattice QCD should take into account the masses of a heavy and a light prop- 
agating quark. The step scaling method faces the challenging requirements 
discussed above by adopting a two-step strategy. As a first step, the decay 
constant is computed on a small volume, where the light quark is squeezed, 
and the heavy quark propagates with a high resolution. At this stage, the 
mass of the heavy quark can be raised up to very large values, with controlled 
discretization errors, and the decay constant can be directly simulated at the 
physical heavy quark mass. As a second step, the finite size effects of this 
calculation are removed by evolving the decay constant toward large volumes. 
The evolution is realized according to the identity 



where the basic ingredient is the ratio of the decay constants computed on 
two different volumes at the same values of the mass parameters 



Throughout the paper we refer to the step scaling function in the continuum 
limit as a (greek lowercase) and to the step scaling function at finite lattice 
spacing as E (greek uppercase). 

This quantity represents a non-perturbative calculation of the finite vol- 
ume effects. In principle, its dependence on the quark masses can be very dif- 
ferent from the one of the decay constants themselves. In effects, it has been 
shown [28] that the cr-ratio's are characterized by a very shght linear depen- 
dence upon the inverse of m/j, due to cancellations of additional heavy quark 
mass dependences between the numerator and the denominator of cq. (2). 
This suggests a concrete way to connect the finite volume decay constant to 
physical volumes: 

• given a couple of physical volumes {Lk-i, Lk) and a finite lattice spacing a, 
the step scaling function is simulated on the lattice for a set of heavy and 
light quark masses. In order to identify the quark mass on a finite volume, a 
RGI quark mass scheme is adopted [30,31] and units are fixed through the 
tq scale [32,33]. Throughout the paper we fix tq = 0.5 fm. The light quark 
masses are kept around the strange mass throughout the whole procedure. 

• A set of different simulations are done at fixed physical volumes {Lk-i, Lk) 
but with different lattice spacings, in order to perform the continuum ex- 



fhe{Loo) = fhi{Lo) 




Lq <. Li <. L2 ■ ■ ■ 



(1) 




(2) 
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trapolation of the step scaling function at given heavy and hght RGI quark 
masses. The ratio s between the two vohimes should be chosen small enough 
to cope with the increase of lattice sites without exceeding computational 
resources. On the other hand, it should be large enough to reach large vol- 
umes in few steps. A value s = 2 is a good compromise. The continuum step 
scaling functions are then linearly extrapolated in the inverse of the RGI 
heavy quark masses up to m^'^^ or m^'-^^ , according to the heavy flavors of 
the meson. 

• As a starting value for the finite volume, we chose to set Lq = 0.4 fm. This 
allows to reach a volume L2 — 1.6 fm, after just two evolution steps, which 
is adequate to accommodate the heavy-light mesons at the physical values 
of the light quark masses. 

In order to match subsequent steps, the knowledge of the bare coupling 
go{ci) as a function of the lattice spacing is required at very small couplings. 
The problem has been recently addressed in [34] and solved by a renormaliza- 
tion group analysis. 

3 Observables 

The step scaling function is calculated within the SF [35,36], which has al- 
ready been apphed to a number of different finite size problems [37,31,38,39,24] . 
The lattice topology is T x with periodic boundary conditions on the space 
directions and Dirichlet boundary conditions along time. We use the following 
set of parameters 

r = 2L, C = C" = 0, ^ = (3) 

where C and C represent the boundary gauge fields and ^ is a topological 
phase which affects the periodicity of the fermion in the space directions. Lat- 
tice discretization is performed using non-perturbative 0{a) improved clover 
action [40] and operators. In order to set the notation, let 

P{x)^i^,{x)j5i^j{x) (4) 

be the axial current and the pseudoscalar density {i and j are fiavor indices) . 
The improvement of the axial current is obtained through the relations 

A'^{x) = A^{x) + acA d^P{x) (5) 

where = {8^ + d^)/2 and 5^, 9* are the usual forward and backward lattice 
derivatives respectively. For what concerns the improvement coefficients ca, 
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we use the non-perturbative results of [40]. The correlation functions used 
to compute the meson decay constants are defined by probing the previous 
operators with appropriate boundary quark sources 



•^iK) = -^E(C,(y)750(z)4(^)) 



•^pW = -^E(C,(y)75C.(z)P(x)) 



(6) 



where Cj(y) diy) can be considered as quark and anti-quark boundary 
states. 

The renormalization of the axial current is realized according to the fol- 
lowing relation 



A^(xo) ^Za(1 + bA am) £^(xo) 



(7) 



here am is the bare quark mass defined as 



am,- 



1 1 

ki kc. 



(8) 



The renormalization constant Za has been computed non perturbatively in 
[41]. For the improvement coefficient Ba we use the perturbative results quoted 
in [42] (at the values of the bare couphng, /3 ~ 7.0, used in the numerical sim- 
ulations the one-loop contribution to 6a differs from the tree-level of 10%). 



The so-called bare current quark masses are defined through the lattice 
version of the PCAC relation 



m, 



wi 

ij 



doJ^Ajxp) + acAdQdoJ^p{xo) 
2J^p(xo) 



(9) 



These masses are connected to the renormalization group invariant (RGI) 
quark masses, according to the definitions given in [30], through a renormal- 
ization factor which has been computed non-perturbatively in [31]: 



mg^^ = ZM{go) 



l + {bA- bp) 



ami + amj 



<'(^o) 



(10) 



where am^ is defined in eq. (8). 

The combination Ba — bp of the improvement coefficients of the axial 
current and pseudoscalar density has been non-perturbatively computed in 
[43,44]. The factor ZM{go) is known with very high precision in a range of 
inverse bare couplings that does not cover all the values of /3 used in our 
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simulations. We have used the results reported in table (6) of ref. [31] to 
parametrize ^m(5'o) in the enlarged range of P values (5.9, 7.6). 

The RGI mass of a given quark is obtained from eq. (10) using the diagonal 
correlations 

- rn^' (11) 

Prom non-diagonal correlations in eq. (10) one obtains different 0{(i) improved 
definitions of the RGI i-quark mass for different choices of the j-fiavor: 

-S' = - (12) 

All these definitions must have the same continuum limit because the depen- 
dence upon the j-fiavor is only a lattice artifact. Further, for each definition 
we use in eq. (9) either standard lattice time derivatives as well as improved 
ones [43,44]. 

Another non-perturbative 0{a) improved definition of the RGI quark 
masses can be obtained starting from the bare quark mass 

mf^^ = ZM{go) Z{go) [1 + 6^ arrii] rrii (13) 
where the improvement coefficient and the renormalization constant 

Ztoo) = %^ (14) 

have been non-perturbatively computed in ref. [43,44]. 

Equations (11), (12) and (13) give us different possibilities to identify the 
valence quarks inside a given meson (fixed by the values of the bare quark 
masses). The procedure is well defined on small volumes because the RGI 
quark mass is a physical quantity that does not depend upon the scale, given 
in the SF scheme by the volume, and is defined in terms of local correla- 
tions that do not suffer finite volume effects. Each pair {mf^^ ^mf^^^ fixed 
a priori is matched, changing the values of the hopping parameters, by the 
different definitions of equations (11), (12) and (13), and leads to values of the 
corresponding decay constants differing by O(a^) lattice artifacts. We take 
advantage of this plethora of definitions by constraining in a single fit the 
continuum extrapolations (see Figure [1,3,6]). 

The meson masses are extracted from the so-called effective mass 

aMxixo) = ^ In [J^x{xo - a)/J^x{xo + a)] (15) 
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where Tx is one of the correlations defined in (6). Correspondingly, the meson 
decay constants are defined as 

JM = / = (16) 

^L^Ma(T/2) V^i 

where is the boundary-to-boundary correlation needed in order to cancel, 
in the ratio, the renormalization constants of the boundary quark fields: 

= E (C,(y)750(z)C.(u)75Ci(w)) (17) 

y,z,u,w 

We want to stress that our choice of defining the decay constant in the middle 
of the lattice, xq = T/2, does not introduces other length scales than L into 
the calculation. Indeed we take T — 2L and the step scaling technique (see 
eq.(l)) connects xq — L^nin, where the decay constant has been defined on the 
smallest volume, with xq — L^^x, where one expects to be free from finite 
volume effects. 



4 Numerical simulations 



In this section we report, step by step, the results of the calculations. 
4-1 Small volume decay constants (Lq — 0.4 im) 

Simulations of the decay constants on the smallest volume {Lq — 0.4 fm) 

have been performed at five different lattice spacings using the geometries 
24 X 12\ 32 X 16■^ 40 x 20'1 48 x 24'"^ and 64 x 32^. For each discretization, a 
set of eight quark masses have been simulated. Two of the heavy masses have 
been chosen around the bottom quark mf^^ = 6.73(16) GeV [29]. Other two 
have been chosen in the region of the charm quark mf^^ = 1.681(36) [29]. An 
additional heavy quark has been simulated with mass 4.00 GeV. Three light 
quark have been simulated with RGI masses of 0.14 GeV, 0.10 GeV and 0.06 
GeV. Using the accurate determination of the RGI strange quark mass given 
in [45] we have fixed one of the simulated fight quarks to be the physical s. We 
will combine this finite volume calculations with the ones of the step scaling 
functions to provide results for the heavy-light decay constants with light 
quarks around the the strange in the continuum and on the large volume. All 
the parameters of the five different simulations are summarized in Table [1]. 

We have obtained different set of data by using the different definitions 
of the RGI quark masses given in the equations (11), (12) and (13). The con- 
tinuum results are thus obtained trough a combined fit of all the set of data. 
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Fig. 1. Continuum extrapolation on the small volume of /^^(Lo). The different values 
of the decay constant correspond to different definitions of the RGI quark masses given 
in equations (11), (12) and (13). For each definition the two points at largest lattice 
spacing has been shown but not included in the fit. Units are in GeV. 

linear in (a/ro)^, as shown in Figure [1] in the case of the bs meson. For 
each set of data we have included in the fit the three points nearest to the 
continuum obtaining a global x^/dof = 0.89 to be compared with the x^s of 
each individual definition listed in the figure. At this small volume, we are 
in a region of small bare couplings (^^q ~ 0.85) where it is legitimate to use 
perturbative values for the improvement coefficient b^- The systematics intro- 
duced in the calculation by the continuum extrapolations have been estimated 
repeating the fits linear in (a/ro)^ including, for each set of data, only the two 
points nearest to the continuum. We find a deviation leading to a correspond- 
ing systematic error of the order of 1% that will be given to the results on 
the large volume added in quadrature with an error of about 2% coming from 
the uncertainties on the lattice spacing and on the renormalization factors. 
The latter have been evaluated by moving the points as a consequence of the 
change, within the errors, of the lattice spacings and of the renormalization 
constants and by repeating the whole analysis. 

The numbers we obtain are: 

fsALo) = ^75{2)MeV JdALq) = QU{3)MeV (18) 

The errors quoted at this stage are statistical only, evaluated by a jackknife 
procedure. 

Due to the compression of the low energy scale, these results are higher than 
the large volume ones obtained after the step scaling functions multiplication 
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Fig. 2. The figure shows the step scaling functions S(Lo) as functions of mf*^^, 



for 



the simulation of the first evolution step corresponding to (3 = 6.963. The different sets 
of data correspond to the values of 7112'^^ . As can be seen the step scaling functions 
approach a plateau for high values of mf^^. Similar plots can be obtained for the other 
values of the bare couplings. 
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Fig. 3. Continuum extrapolation on the first evolution step, Lq Li, of the step 
scaling function, S(Lo), of the pseudoscalar meson corresponding to the heavy quark 
of mass mi'^^ = 3.55 GeV and to the s-quark. The two sets of data are obtained using 
the two definitions of RGI quark masses of equations (11) and (13). Units are in GeV. 
Similar plots can be obtained for the other combinations of quark masses used in our 
simulations. 



chain (see eqs. (21), (22) and (23)). 



4.2 First evolution step (Lq — > Li) 

The finite volume effects on the decay constants calculated on Lq, are 
measured by doubling the volume, Li = 0.8 fm, and by using the step scaling 
function of eq (2). 

The continuum extrapolations have been obtained by simulating the step 
scaling functions with three different discretizations of Lq, i.e 16 x 8'^, 24 x 12^ 
and 32 x 16'^. The volume Li has been simulated starting from the discretiza- 
tions of Lq, fixing the value of the bare coupling and doubling the number of 
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Fig. 4. The figure shows the continuum extrapolated step scaling functions a{Lo) as 
functions of . The heavy extrapolations are shown only for the hs set of data. 

lattice points in each direction. 

The simulated quark masses have been halved with respect to the masses 
simulated on the small volume in order to have the same amount of the dis- 
cretization effects proportional to am. The set of parameters for the simula- 
tions of this evolution step is reported in Table [2] . 

The step scaling functions at f3 = 6.963 are plotted, at fixed m^*^^, as func- 
tions of mf"^-^ in Figure [2]. As can be seen, S(Lo) is almost flat in a region of 
heavy quark masses starting around the charm mass. The hypothesis of low 
sensitivity upon the high-energy scale is thus verified. The value of the step 
scaling functions for the s quark are obtained trough linear interpolation. 
In Figure [3] are reported the results of the continuum extrapolation of the 
step scaling function, T,{Lq), of the pscudoscalar hs meson corresponding to 
the heaviest quark simulated in this step (m^'^^ = 3.55 GeV). 
The residual heavy mass dependence of the continuum extrapolated step scal- 
ing functions is very mild, as shown in Figure [4] in the plot of as, as a function 
of the inverse quark mass. The continuum results are linearly extrapolated at 
the values of the heavy quark masses used in the small volume simulations. 

The numbers we get at this step are: 

(Lo) = 0.417(3) azj, (Lq) = 0.414(3) (19) 

The step scaling functions are free from the systematic errors coming from 
uncertainties on and since the multiplicative improvement and renor- 
malization factors cancel exactly in the ratio, being the numerator and the 
denominator evaluated at the same lattice spacing. 

4-3 Second evolution step (Li — >• L2) 

In order to have the results on a physical volume, L2 = 1.6 fm, a second 
evolution step is necessary. This is done computing the step scaling functions 
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Fig. 5. The figure shows the step scaling functions as functions of mf , for 

the simulation of the first evolution step corresponding to /? = 6.420. The different sets 
of data correspond to the values of 7112'^^ . As can be seen the step scaling functions 
approach a plateau for high values of mf^^. Similar plots can be obtained for the other 
values of the bare couplings. 
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Fig. 6. Continuum extrapolation on the first evolution step, Li 1-^ L2, of the step 
scaling function, Ti{Li), of the pseudoscalar meson corresponding to the heavy quark 
of mass mf^^ = 2.00 GeV and to the s-quark. The two sets of data are obtained using 
the two definitions of RGI quark masses of equations (11) and (13). Units are in GeV. 
Similar plots can be obtained for the other combinations of quark masses used in our 
simulations. 

of eq. (2) at L = Li, by the procedure outlined in the previous section. The 
parameters of the simulations are given in Table [3]. 

Also in this case, the values of the simulated quark masses have been 
halved with respect to the previous step, owing to the lower values of the sim- 
ulation cutoffs. Even if we are lowering again the values of the quark masses, 
the linear extrapolations at the values of the heavy quark masses used on the 
small volume appears to be still valid and under control; see Figure [5,7]. 
The value of the step scaling functions for the s quark are obtained trough 
linear interpolation. 

Figure [6] shows the continuum extrapolation of the step scaling function, 
S(Li), of the hs meson corresponding to the heaviest quark simulated in this 
step {m^^^ = 2.00 GeV). 



11 




1/m! 



RGI 



Fig. 7. The figure shows the continuum extrapolated step scaling functions (t{Li) as 
functions of l/mf^^. The heavy extrapolations are shown only for the hs set of data. 
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Fig. 8. Chiral extrapolation of the continuum results for the pseudoscalar decay con- 
stants on the large volume. Units are in GeV. 



The numbers for this step are: 

aB, (Li) = 0.97(3) (L,) = 0.90(2) 



(20) 



5 Physical results 

In this section we combine the results of the small volume with the re- 
sults of the step scaling functions to obtain, according to eq. (1), the physical 
numbers. In the end, we get: 

/b. = 192(6)(4) MeV fo, = 240(5)(5) MeV (21) 

The first error is statistical while the second one is our estimate of the system- 
atics due to the uncertainties on the continuum extrapolations, on the scale 
and on the renormalization factors, as already discussed in sec. 4.1. 
Note that our value for fo^ agrees with the average of dedicated calculations 
performed on large volumes [15,16]. This validates our choice of stopping at 
L = 1.6 fm that, of course, can be explicitly checked to be safe by performing 
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further evolution steps. 

Using the strategy outhned in the previous sections, we have calculated also 
the decay constant of the Be meson. The result we obtain is 

/b, = 347(5) (8) MeV (22) 



that represents the first determination of this quantity from quenched lattice 
QCD. 

The chiral behavior of heavy-hght pseudoscalar decay constants has been 
shown [46,47,48,49] to contain logarithmic terms (x-logs) that are diverging 
in the chiral limit, at variance with the unquenched case where these terms 
only affect the form of the extrapolation. 

In Figure [8] we show the chiral extrapolations for the continuum heavy-light 
pseudoscalar decay constants. The data corresponding to the different values of 
the parameter A have been extrapolated using the parametrization suggested 
in [49]. As can be seen from the figure, the presence of the unphysical quenched 
X-logs make the extrapolations unreliable down to the li-quark mass while the 
strange region seems to be dominated by a linear behavior. Nevertheless, in 
the literature values extrapolated linearly in the light quark mass have been 
quoted. 

For a historical comparison we can quote our own: 



j^Hnear ^ 171(8) (4) MeV /j^"^"'' = 221 (7) (5) MeV 

= 1-12(2)(1) 7^ = 1-09(1)(1) (23) 

JB Jd 

that differ by a large factor from the unreliable values obtained from the fits 
shown in Figure [8] because of the diverging x-logs. The numbers quoted 
above should then be retained for an historical comparison only and should 
not be quoted as the results of the quenched approximation. 



6 Conclusions 



In this work we have calculated the i^^-meson decay constant in the con- 
tinuum hmit of quenched lattice QCD. The results are obtained through a 
finite volume recursive procedure where the heavy quark masses have been 
obtained, using the same method, in a previous work. 

The main achievement of this computation with respect to our previous 

determination of the same quantities are the extrapolations to the continuum 
limit. Our systematic errors are due to the extrapolations to the continuum 
and to the physical heavy quark masses. An additional unknown systematics 
comes from the quenched approximation that is believed to produce visible 
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effects on the meson decay constants. The known pathological behavior of 
quenched x-QCD does not allow us to quote a quenched value for fs- 

Upcoming new powerful super-computers will make affordable straight 
calculations of lattice 6-physics without recursive methods but still in the 
quenched approximation. In this scenario our method could provide the oppor- 
tunity of studying txnquenched 6-physics and/or to deal with other demanding 
two-scales problems. 
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Lo/a fee 


fe 


m««5^ (GeV) 






0.115528 


7.51(9) 






0.116762 


6.91(8) 






0.123555 


4.019(43) 


6.737 


12 0.13520(1) 


0.130384 


1.702(20) 






0.130089 


1.604(18) 






0.134801 


0.1347(31) 






0.134925 


0.0929(29) 






0.135048 


0.0513(29) 






0.120081 


7.14(8) 






0.120988 


6.63(7) 






0.126050 


4.024(44) 


6.963 


16 0.134827(6) 


0.131082 


1.696(19) 






0.131314 


1.591(18) 






0.134526 


0.1381(30) 






0.134614 


0.0978(28) 






0.134702 


0.0574(28) 






0.122666 


7.03(11) 






0.123437 


6.53(10) 






0.127605 


3.97(6) 


7.151 


20 0.134492(5) 


0.131511 


1.716(27) 






0.131686 


1.617(25) 






0.134277 


0.1257(36) 






0.134350 


0.0829(33) 






().l:-!il22 


().()l()7(:-!2) 






0.124176 


7.11(8) 






0.124844 


6.61(20) 






0.128440 


4.018(44) 


7.300 


24 0.134235(3) 


0.131800 


1.695(19) 






0.131950 


1.592(18) 






0.134041 


0.1374(27) 






0.134098 


0.0971(24) 






0.134155 


0.0567(24) 






0.126352 


7.10(8) 






0.126866 


6.60(7) 






0.129585 


4.016(44) 


7.548 


32 0.133838(2) 


0.132053 


1.698(19) 






0.132162 


1.595(18) 






0.133690 


0.1422(27) 






0.133732 


0.1021(25) 






0.133773 


0.0618(23) 



Table 1 

Simulation parameters at Lq = 0.4 fm. The RGI quark masses are obtained using 
eq. (11). 
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Lo/a fee 


fe 


m'^'^' (GeV) 






0.120674 


3.543(39) 






0.122220 


3.114(34) 






0.126937 


1.927(21) 


6.420 


8 0.135703(9) 


0.134304 


0.3007(36) 






0.134770 


0.2003(28) 






0.135221 


0.1028(21) 






0.1249 


3.542(39) 






0.1260 


3.136(34) 






0.1293 


1.979(22) 


6.737 


12 0.135235(5) 


0.1343 


0.3127(38) 






0.1346 


0.2090(28) 






0.1349 


0.1080(21) 






0.127074 


3.549(39) 






0.127913 


3.153(35) 






0.130409 


2.003(22) 


6.963 


16 0.134832(4) 


0.134145 


0.3134(38) 






0.134369 


0.2112(28) 






O.l:-! 159:-! 


().ll)8(){20) 



Tabic 2 

Simulation parameters for the first evolution step Lq — > Li = 0.8 fm. The RGI quark 
masses are obtained using eq. (11). 
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Li/a fee 


fe 


m'^'^' (GeV) 






0.118128 


2.012(22) 






0.121012 


1.551(17) 






0.122513 


1.337(15) 


5.960 


8 0.13490(4) 


0.131457 


0.3154(36) 






0.132335 


0.2322(28) 






0.133226 


0.1466(44) 






0.124090 


1.984(22) 






0.126198 


1.584(17) 






0.127280 


1.389(15) 


6.211 


12 0.135831(8) 


0.133574 


0.3493(39) 






0.134177 


0.2550(29) 






0.134786 


0.1510(19) 






0.126996 


1.933(21) 






0.128646 


1.547(17) 






0.129487 


1.355(14) 


6.420 


16 0.135734(5) 


0.134318 


0.3016(34) 






0.134775 


0.2038(24) 






().i;r.23r> 


0. 1055(15) 



Tabic 3 

Simulation parameters for the first evolution step Li — > L2 = 1.6 fm. The RGI quark 
masses are obtained using eq. (11). 
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